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1.  Introduction 

This  report  updates  the  radiation/energy  balance  model.  [1]  Most  of  the  original 
model  concepts  and  equations  are  retained;  however,  a  few  significant 
modifications  and  additions  have  been  made.  The  changes  include  an  expanded 
and  corrected  table  of  soil  thermal  constants,  tables  of  surface  roughness,  albedo 
values,  [2,3],  and  a  dramatically  improved  method  for  bridging  the  fluxes  through 
the  sunrise  neutral  period.  Finally,  we  have  corrected  editorial  errors  in  the 
earlier  report. 

The  optical  turbulence  structure  parameter,  C.,  appears  in  equations  used  to 
estimate  the  effects  of  atmospheric  turbulence  on  imagery  and  electromagnetic 
(EM)  propagation.  For  many  optical  systems  corresponds  with  degradation 
of  performance.  [4,5,6]  Basic  equations  of  [5,7,8,9,10,11]  include  the  real 
index  of  refraction  gradient  as  a  coefficient,  which  is  a  function  of  the 
temperature  and  moisture  gradients  (preferably  potential  temperature  and  specific 
humidity  according  to  Tatarski  [5]).  Temperature  and  ihoisture  gradients  can  be 
approximated  fi'om  sensible  and  latent  heat  flux  estimates,  and  the  fluxes  can  be 
obtained  from  radiation/energy  balance  formations. 

Numerous  radiation/energy  balance  models  appear  in  the  literature 
[12,13,14,15,16,17,18]  varying  from  comparatively  simple  to  academically 
complex  and  requiring  different  amounts  and  numbers  of  inputs  and  computer 
capabilities.  This  report  presents  a  model  (hereafter  known  as  the  RT  model)  that 
was  developed  for  imaging  and  EM  propagation  applications  when  minimum 
atmospheric  information  is  available.  This  report  also  provides  the  concept  and 
equations  that  make  up  the  second  formali^  version  of  the  model.  Other 
reports  will  follow  as  the  model  evolves. 

The  primary  ouq)uts  of  the  RT  model  are  estimates  of  sensible  and  latent  heat 
fluxes,  which  in  turn,  are  used  for  estimating  gradients  of  potential  temperature 
and  specific  humidity,  making  use  of  Monin-Obukhov  similarity  relations.  These 
gradients  are  then  used  to  estimate  the  optical  turbulence  parameter  needed  in 
the  imaging  and  propagation  equations. 

The  remainder  of  this  report  is  structured  into  the  following  sections:  Model 
Concq}t,  Model  Equations,  Computational  Procedure,  Model  Results,  Discussion, 
Summary  and  Conclusions,  and  References. 
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2.  Model  Concept 


The  following  model  concept  is  takra  from  the  operational  scenario.  An  interest 
is  developed  in  estimating  the  optical  turbulence  at  a  site  for  different  times 
of  a  day.  The  day  of  interest  is  known,  and  the  longitude  and  latitude  of  the  site 
are  known.  The  soil  type  may  be  determined  from  the  USDA  Soil  Conservation 
Service’s  soil  maps.  Soil  wetness  (dry,  moist,  or  wet)  may  be  determined  from 
meteorological  reports.  However,  if  the  soils  taxonomy  and  meteorological 
rqwrts  are  not  available,  on-the-spot  judgments  are  made  by  examining  samples 
of  soil  from  the  surface  to  10  cm  below  the  surface.  A  judgment  of  sky 
conditions  is  also  made;  in  particular,  amount  of  cloud  cover  (in  tenths  Ax)m  0 
to  1),  cloud  height,  and  density  or  sun  blocking  capacity  of  the  clouds  (0  to  3). 
From  the  information,  estimations  are  made  for  net  radiative  flux,  sensible  heat, 
ground  heat,  and  latent  heat  fluxes  using  the  convention  and  equations  of  radiation 
and  energy  balance  proposed  by  Carson.  [18] 


2.1  Radiation  Balance  Concepts 


Carson  [18]  states  that  the  net  radiation  flux  at  the  soil  surface  is  equal  to  the 
sum  of  the  net  short-wave  radiative  flux,  Rsn,  and  the  net  long-wave  radiative 
flux,  Ru,; 


^Sff  *  ^Ut' 


The  short-  and  long-wave  fluxes  are  illustrated  in  figure  1. 


SHORT-WIVVE  FLUXES 


LONG-WWfE  FLUXES 
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Figure  1.  Schematic  representation  of  short-wave  and  long- wave  radiative  fluxes  at  a 
hare-soil  surface. 
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In  figure  1,  R,t  is  ihc  downward  short-wave  radiative  flux  (including  both  the 
direct  solar  flux  and  diffuse  radiation  from  the  sky).  However,  some  of  the 
short-wave  radiation  is  reflected  at  the  surface  of  the  earth  so  that 

^sN  “  ■  “  ^s*  * 


where 

a  =  the  surface  short-wave  reflectivity  (albedo). 

Albedo  is  not  as  simple  as  it  appears  in  equation  2;  it  is  a  function  of  soil  type, 
color,  vegetative  cover  and  the  elevation  angle  of  the  sun.  Values  of  a  as  a 
function  of  soil  and  vegetation  are  given  in  tables  la  through  Ic,  mainly  from  an 
extensive  compilation  by  Hansen  [2]  and  Pielke  [IS];  the  functional  form  of  a, 
relative  to  the  position  of  the  sun,  is  given  in  section  3  (Model  Equations). 

Carson  [18]  notes  that  if  is  the  downward  long-wave  radiation  and  a  is  the 
surface  absorptivity  to  long-wave  radiation,  the  net  incoming  flux  from  the 
atmosphere  is  (a  Rl^).  Using  Stefan’s  law,  the  upward  flux  caused  by  thermal 
emission  at  the  surface  of  the  earth  can  be  written  as 

c  <rTj 


where 

T,  =  the  effective  surface  temperature 

6  =  the  long-wave  emissivity  at  the  surface 

a  =  the  Stefan-Boltzman  constant. 

Hence,  the  net  long-wave  radiative  flux,  Rln,  is 

(3) 

where 

a  =  5.6697  x  10*. 

It  is  common  practice  to  combine  the  definition  of  e  with  Kirchhoffs  law  so  that 
a  =  €.  However,  c  varies  with  soil  type,  vegetation,  and  snow  or  water  cover 
as  given  in  table  2  from  Pielke.  [IS] 


Table  la.  Typical  values  of  albedo,  a,  for  various  surfaces 


Surfsu:e  a 


Snow  and  Ice 

Snow,  fresh  fallen  75-95 

Snow,  thawing  30-65 

Snow,  old  40-70 

Snow,  icy  75 

Ice,  grey  60 

Ice,  white  75 

Ice,  water  covered  26 

Ice,  light  snow  cover  31 

Ice,  porous  and  melting  41 

Desert  shrublands,  snow  covered  18-19 

Conifer  forest,  snow  covered  59-67 

Mixed  forest,  50  cm  snow  covered  20 

Grasslands,  snow  covered  46-50 

Crops,  snow  covered  18-19 

Tundra,  snow  covered  59-67 

Soils  and  Rocks 

Soil,  dark,  plowed,  wet  6 

Soil,  dark,  plowed,  dry  8 

soil,  light,  plowed,  wet  8 

Soil,  light,  plowed,  dry  16 

Soil,  dark,  wet  8 

Soil,  dark,  dry  13 

Soil,  light,  wet  10 

Soil,  light,  dry  18 

Dark  organic  soils  10 

Dark  grey  silt  12 

Red  soils  17 

Brown  soils  17 

Clay,  wet  16 

Clay,  dry  23 

Lo^,  wet  16 

Loam,  dry  23 


Table  lb.  Typical  values  of  albedo,  a,  for  various  surfaces 
(continued) 


Suiface 


a 


Soils  and  Rocks  (continued) 


Gay  loam,  wet 

19 

Gay  loam,  dry 

13 

Sandy  soil,  wet 

20 

Sandy  soil,  dry 

25 

Sand,  white,  wet 

25 

Sand,  white,  dry 

35 

Peat  soils 

5-15 

Lime 

45 

Gypsum 

55 

Lava 

10 

Granite 

12-18 

Rock,  wet 

20 

Rock,  dry 

35 

Stone 

30 

Urban  (People-Influenced) 

Road,  blacktop 

14 

Road,  stone 

15 

Road,  dirt,  wet 

18 

Road,  dirt,  dry 

35 

Road,  clay,  wet 

20 

Road,  clay,  diy 

30 

Road,  asphalt,  wet 

10 

Road,  asphalt,  dry 

15 

Parking  lot,  black  top 

8 

G>ncrete,  new,  white 

37 

Buildings 

9 

Developed  urban  area  (average) 

15 

Roof,  thatched,  new 

20 

Roof,  thatched,  old 

15 

Roof,  tiled,  dirty 

8 

Table  Ic.  Typical  values  of  albedo,  a,  for  various  surfaces  (continued) 


Suiface  a 


Crops,  Natural  Terrain  and  Vegetation 

Fallow  field,  wet  5-7 

Fallow  field,  dry  8-12 

Spring  wheat  10-25 

Winter  wheat  16-23 

Rice  paddy  12 

Sugar  cane  15 

Cocoa  16 

Ground  nuts  17 

Winter  rye  18-23 

Beets  18 

Maize  18 

Tt^acco  19 

Potatoes,  Yams  19 

Alfalfa  23-32 

Cotton  20-22 

Sorghum  20 

Lettuce  22 

Forest,  coniferous  5-15 

Forest,  deciduous  10-20 

Grass,  green  26 

Meadows,  green  10-20 

Coniferous  trees,  dormant  12 

Deciduous  trees,  dormant  12 

Tall  grass,  dormant  13 

Mowed  grass,  dormant  19 

Tundra  15-20 

Savanna  15 

Steppe  20 

Sand  dune,  wet  20-30 

Sand  dune,  dry  35-45 


Table  2.  Typical  values  of  emissivity,  c,  for  various  surfaces 


Ground  Cover 

e 

Fresh  snow 

0.99 

Old  snow 

0.82 

Dry  sand 

0.95 

Wet  sand 

0.98 

Dry  peat 

0.97 

Wet  peat 

0.98 

SoUs 

0.9-0.98 

Asphalt 

0.95 

Concrete 

0.71-0.9 

Tar  and  gravel 

0.92 

Limestone  gravel 

0.92 

light  sandstone  rock 

0.98 

Desert 

0.84-0.91 

Grass  lawn 

0.97 

Grass 

0.90-0.95 

Deciduous  forests 

0.95 

Coniferous  forests 

0.97 

Urban  area  (range) 

0.85-0.95 

2.2  Surface  Energy  Flux  Balance  Concepts 

Carson  [18]  writes  the  energy  flux  balance  at  the  soil  surface  as 

-  H  *  fE  ♦  G,  (4) 

where 

Rn  =  the  net  radiative  flux 

H  =  the  turbulent  sensible-heat  flux 

L’E  =  the  latent-heat  flux  caused  by  surface  evaporation 

G  =  the  flux  of  heat  into  the  soil. 

The  sign  (directional)  convention  of  equation  4  is  shown  in  figure  2. 

Equations  for  H,  L’E,  and  G  are  given  in  section  3  (Model  Equations).  Having 
established  the  basic  radiation/energy  balance  equations  and  their  directional 
conventions,  the  equations  chosen  to  evaluate  them  are  discussed.  Although  a 
multitude  of  relations  exist  to  choose  from,  the  expressions  that  best  satisfy 
constraints  necessary  for  this  application  are  used. 


ENERGY  FLUX  BALANCE 


3;  Model  Equations 

3.1  Short-Wave  Solar  Radiation 

The  equations  used  to  compute  the  incoming  short-wave  solar  radiation  for 
cloudless  skies  are  pattern^  after  Meyers  and  Dale  [19]  and  Miller  and 
Ricklin.  [4]  The  equations  are  augmented  with  empirical  results  by  Haurwitz  [20] 
to  account  for  cloudy  skies. 

For  clear  skies  Meyers  and  Dale  [19]  write 

Rg^  =  7  =  /„  Tg  cos  Z,  (5) 

where 

I,  =  the  extraterrestrial  flux  density  at  the  top  of  the  atmosphere  on  a 
surface  normal  to  the  incident  radiation 
Z  =  the  solar  zenith  angle 

Tk  =  the  transmission  coefficients  for  Rayleigh  scattering 
T,  =  absorption  by  permanent  gases 

T*  =  water  vapor 

T,  =  absoiption  and  scattering  of  aerosols. 

The  incident  radiation,  ]«,  changes  throughout  the  year  because  of  changes  in  the 
earth-sun  distance  and  is  adjusted  by  using  the  equation 

/„  =  1353  (W  m-2)  {  1  +  0.034  cos  [  2t  (n’  -  1  )  /  365  ]  }, 


where 

n’  =  the  Julian  day. 

The  solar  zenith  angle  is  computed  by  using 

Z  =  cos’*  {  sin  (S)  sin  (D)  +  cos  (6)  cos  (D)  cos  (77’)}, 

where 

5  =  the  latitude 

D  =  the  declination  angle 

H’  =  the  solar  hour  angle. 

Miller  and  Ricklin  [4]  compute  the  solar  declination  angle  using  equations  by 
Woolf  [21]: 


15 


sin  Z>  B  sin  (23.4438)  sin 


(8) 


where  in  degrees  is 

/8  =  7  +  0.4087  sin  (7)  +  1.8724  cos  (7)  -  0.0182  sin  (27) 
0.0083  cos  (2y) 

and  7  in  degrees  is 

7  =  279.9348  +  d. 


(9) 

(10) 


The  angle  d  is  the  angular  fraction  of  a  year  represented  by  a  particular  date  and 
may  be  calculated  by 

d  -  {number  of  day  of  year  -  1)  (360/365.242). 


The  solar  hour  angle  H’  in  degrees  is  a  measure  of  the  longitudinal  distance  from 
the  sun  to  the  point  of  calculations  given  by 

IT  15  (T  -  if)  -  11, 


where 

T  =  the  Greenwich  mean  time  (GMT)  of  the  calculations 
M  =  the  time  in  hours  after  midnight  of  the  passage  of  the  sun  over  the 
Greenwich  meridian  or  true  solar  noon 
11  =  longitude,  counted  positive  west  of  Greenwich. 

In  terms  of  d  defined  in  equation  1 1, 

M  *  12.0  +  0.12357  sin  (d)  -  0.004289  cos  (d) 

*  0.153809  sin  (2d)  *  0.060783  cos  (2d) 


The  solar  hour  angle  H’  that  relates  to  sunrise  and  sunset  is  found  using 
equation  7  as 

sin  i4  ~  sin  5  sin  D 
cos  5  cos 

where 

A  =  the  sunrise/sunset  solar  elevation  angle. 

At  ground  level  Woolf  [21]  sets  A  =  -0.9  degree.  The  solar  day  extends  from 
M  -  H’  (sunrise)  to  M  +  H’  (sunset). 


An  empirical  equation  for  T,  in  equation  5  by  Kondartyev  [22]  and  modified 
by  Atwater  and  Brown  [23]  to  include  forward  scattering  is 

7,  7,  «  1.021  -  0.084  [m{9A9  Px  10  *  +  0.051)^,  <^5) 

where 

P  =  the  surface  pressure  (kPa) 

m  s  the  optical  air  mass  at  a  pressure  of  101.3  kPa,  given  by 

m  =  35  (1224  cos^  {Z)  +  l)'*'^ 

where 

Z  the  zenith  angle. 

An  expression  for  computing  the  broad-band  transmission  of  water  vsqwr 
absorption  by  McDonald  [24]  is 

7^  =  1  -  0.077  {wnf  \  (17) 


where 

m  =  the  optical  air  mass  (defined  above). 


The  precipitable  water  vsqior,  u,  is  determined  using  an  expression  by  Smith  [25]; 

u  =  (18) 


«(X  D’ 


where 

P,  =  the  pressure  at  the  surface  of  the  earth 
W,  =  the  mixing  ratio 
g  =  the  acceleration  caused  by  gravity 

and  X  is  given  in  table  3. 

A  simple  treatment  used  by  Meyers  and  Dale  [19]  for  estimating  T,  was  proposed 
by  Houghton  [26] 

7.  =  X",  (19) 


where 

m  =  the  optical  mass 

X  =  an  empirically  derived  constant  (on  the  order  of  0.935). 


\ 
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Table  3.  Seasonal  and  latitudinal  mean  values  of  X 


Season 

Latitudinal  Annual 

Zone  (“  N)  Winter  Spring  Summer  Fall  Average 


0-10 

3.37 

2.85 

2.80 

2.64 

2.91 

10-20 

2.99 

3.02 

2.70 

2.93 

2.91 

20-30 

3.60 

3.00 

2.98 

2.93 

3.12 

30-40 

3.04 

3.11 

2.92 

2.94 

3.00 

40-50 

2.70 

2.95 

2.77 

2.71 

2.78 

50-60 

2.52 

3.07 

2.67 

2.93 

2.79 

60-70 

1.76 

2.69 

2.61 

2.61 

2.41 

70-80 

1.60 

1.67 

2.24 

2.63 

2.03 

80-90 

1.11 

1.44 

1.94 

2.02 

1.62 

Northern 

Hemisphere 

Average 

2.52 

2.64 

2.62 

2.70 

2.61 

To  account  for  the  effect  of  clouds,  a  transmission  coefficient,  derived  empirically 
by  Haurwitz,  [20]  is  introduced.  Haurwitz  computed  the  ratio  of  insolation  widi 
partly  or  completely  covered  sky  to  insolation  of  cloudless  skies,  as  shown  in 
table  4.  Table  4  is  used  to  correct  estimates  of  insolation  for  cloud  cover. 

At  this  point,  equations  are  available  for  computing  the  incoming  short-wave  solar 
radiation  R  which  is  adjusted  for  albedo.  As  indicated  in  table  1,  albedo  is  a 
function  of  soil  type,  color,  moisture,  vegetation,  and  solar  elevation  angle. 
Although  the  first  four  effects  are  reflected  in  table  1,  the  last  is  computed  by 
using  an  approximation  by  Paltridge  and  Platt  [27]; 

a(Z)  =«,•*•(!-  a,)  exp  [-k'  (90“  -  Z)  ],  (20) 


where 

Oi  =  the  small  zenith  angle  value  of  albedo 

Z  =  the  solar  zenith  angle 

and  k’  is  of  order  0.1. 
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3.2  Downward  Long-Wave  Radiation 

Gates  [28]  offers  an  expression  for  downward  long-wave  radiation  for  clear  skies 
based  on  an  empirical  relationship  by  Swinbank.  [29] 

=  -170.9  1.195  a  7^,  (21) 


where 

R  =  watts  per  metei^ 

a  =  5.6697  x  lO* 

T,  =  the  reference  level  temperature  (K). 

Paltridge  and  Platt  [27]  suggest  an  addition  to  equation  21  to  account  for  clouds, 
giving  a  total  expression  of 

»  -170.9  +  1.195  <T  7^  +  0.3  o  t  (cc),  (22) 

where 

Ce  =  the  emissivity  of  the  cloud  base  (given  in  table  5) 

Tc  =  the  temperature  of  the  cloud  base  (K). 

Having  an  estimate  of  the  cloud  height,  can  be  approximated  assuming  an 
average  of  the  dry  and  moist  adiabatic  lapse  rate.  Recent  experience  has  shown 
that  the  empirically  based  constant  in  equations  21  and  22  can  be  adjusted  to 
obtain  improved  radiation  estimates  before  sunrise  and  after  sunset.  Modifying 
the  value  of  the  empirically  based  constant  from  -170.9  to  -130.0  has  been 
somewhat  successful. 


Table  5.  Mean  short-wave  planetary  albedo,  »  and  long-wave  flux  emissivity, 


Cloud  Level 

Cloud  Type 

®ci 

Cci 

1 

Cirrus 

0.35 

0.3 

2 

Alatocumulus-altostratus 

0.55 

0.9 

3 

Low  cloud  2 

0.60 

1.0 

4 

Low  cloud  1 

0.50 

1.0 
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3.3  Upward  Long-Wai^e  Radiation 

The  upward  long-wave  radiative  flux  Ri.t  is  computed  using  an  equation  from 
Yamada  [17]  and  is  expressed  as 

where 

^  -  the  surface  enussivity  (see  table  2) 

(T  -  the  Stefan-Boltzman  constant. 

An  estimate  of  the  effective  ground  temperature  T,  is  required  to  evaluate 
equation  23.  Obukhov  pO]  similarity  theory  parameterizations  are  applied  to 
reach  the  estimate.  For  daytime,  unstable  surhu»  layer  conditions,  it  is  assumed 
that  a  senuempirical  expression  of  the  sensible  heat  flux,  H,  given  by  Angus- 
Lqrpan  [31]  (based  on  the  Penman  p2])  combination  form)  is  adequate  as  a  first 
s^roximation  for  clear  skies  and  dry  ground.  Angus-Leppan’s  expression  is 

H  =  450  C  W  sin  4>„ 

whm 

C  =  the  correction  factor  (figure  3)  for  amount  of  cloud  covn- 

W  =  the  correction  factor  (figure  3)  for  amount  of  ground  wetness 

=  the  solar  elevation  angle  (0,  ~  90  -  Z). 

In  equation  24,  C  =  W  =  1  for  clear  skies  and  dry  ground. 

For  stable  conditions,  night  or  day,  it  is  assumed  that  H  is  constant  and  negative 
in  rign.  Although  this  is  a  relatively  weak  attempt  at  modeling  the  complexities 
of  the  sensible  heat  flux  and  its  associated  surface  layer  temperature  gradient 
structure  during  stable  conditions,  it  allows  the  application  of  a  modding  effort 
to  a  24-h  data  set  with  some  degree  of  confidence.  Literature  suggests  that  in 
many  cases  the  sensible  heat  flux  remains  relatively  flat  or  uniform  in  the  stable 
surface  layer.  [31] 

Having  a  first  estimate  for  H,  the  similarity  forms  for  the  Obukhov  length,  L,  and 
friction  velocity,  u*,  are  computed  using 

C  p  0 

L  =  ^  ^  (25) 

kgH 
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Figure  3.  Estimation  of  parameters  C  and  W. 


where 

Cp  =  the  specific  heat  of  air 

p  =  the  density  of  air 

d„  =  the  reference  level  potential  temperature 

k  =  Karman’s  constant  (0.4) 

g  =  acceleration  caused  by  gravity,  and 

»•  rjui  l^j  »2lan-^|  IL-'j", 

where 

V,  =s  the  reference  level  wind  speed 
j:  =  (1  -  5  i. for  unstable  conditions 

l' 

X  =  1  +  5  4.  for  stable  conditions. 


(26a) 
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For  neutral  conditions 


u*  =  V^k 


(26b) 


In  equations  26a  and  26b, 

Zo  =  the  roughness  length  for  momentum 
d  =  the  displacement  height. 


Approximate  values  for  these  parameters  can  be  estimated  by  Zq  =  0. 14  z,  and 
d  =  0.7  Ze ,  where  z^  is  referred  to  as  the  height  of  the  roughness  element  or  the 
average  height  of  the  terrain  or  ground  cover  (surface  vegetation).  Table  6 
rq)resents  an  extensive  compilation  by  Hansen  [3]  of  typical  roughness  length 
values  for  various  surfaces. 

Iteratively  solving  equations  25  and  26a  or  26b  obtains  values  for  L  and  u*  for 
the  approximated  value  of  H.  From  similarity  forms,  it  can  be  written 

e-  =  (27) 

kgL 


Hence,  a  first  estimate  for  0*,  the  potential  temperature  scaling  constant,  can  be 
computed  and  used  to  approximate  T,  because 


T,  =  m  - 


(28a) 


where 

y  ®  (1  -  15  y)*'^  for  unstable  conditions 
y  =  1  +  5  y  for  stable  conditions. 


For  neutral  conditions 


(28b) 
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Table  6a.  Typical  values  of  roughness  length,  Z„  for  various  surfaces 


Type  of  Surface  Zq  (cm) 


Ice  0.001 

Smooth  mud  flats  0.001 

Dry  lake  bed  0.003 

Tundra,  snow  covered  0.01 

Tundra,  patchy  snow  0.03 

Tundra,  after  snow  melt  0.40 

Tundra,  mid-summer  2.40 

Snow  cover  (Antarctic)  0.01 

Calm  open  sea  0.01 

Desert,  smooth  0.03 

Grass,  closely  mowed  0. 10 

Grass,  short  0.14 

Farmland,  snow  covered  0.20 

Bare  soil,  tilled  0.20-0.60 

Nebraska  prairie  0.70 

Sparse  grass,  10  cm  high  0.70 

Thick  grass,  5-6  cm  high  0.75 

Grassy  plains,  level  topography  1.0 

Kansas  prairie  1.0 

Low  shrubs,  level  topography  2.6 

Grasslands,  18  cm  high  2.7 

Uncut  grass  w/isolated  trees  3.0 

Grass  and  trees,  mixed  3.5 

Sparse  brush,  semi-arid  5.0 

Sparse  grass,  50  cm  high  5.0 

Thick  grass,  50  cm  high  9.0 

Thick  grass,  60-70  cm  high  8.10-15.0 

Brush,  scrub  growth,  open  16.0 

Brush,  scrub  growth,  dense  25.0 

Wooded  country,  level  topography  40.0 

Forested  plateau,  level  topography  70.0-120.0 

Forested  plateau,  rolling  hills  120.0-130.0 

Low  mountains,  hills,  unforested  75.0 

Subtropical  Savannah,  grass,  few  trees  31.0-41.0 

Subtropical  Savannah,  grass  w/many  shrubs  51.0-61.0 

Coniferous  forest  1 10.0 
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Table  6b.  Typical  values  of  roughness  length,  Z,,  for  various  surfaces 
(continued) 


Type  of  Surface 

Zo  (cm) 

Alfalfa 

2.7 

Cashew  orchard,  2  m  high 

3.5-4.0 

Potatoes,  60  cm  high 

4.0 

Farmland,  few  trees 

6.0 

Farmland,  many  hedges 

8.0 

Bean  crop,  1.2  m  high 

7.4 

Cotton,  1.27  m  high 

13.0 

Wheat 

22.0 

Citrus  orchard 

31.0-40.0 

Com,  2.2  m  high 

74.0 

Blacktop,  concrete 

0.002 

Airport  runways 

3.0 

Highways,  railways 

50.0 

Villages,  towns 

40.0-50.0 

Residential,  low  density 

110.0 

City  park 

130.0 

Urban  buildings,  business  district 

175.0-320.0 

High  rise  sqxutments 

370.0 

However,  to  evaluate  equation  28a  or  28b  an  estimate  of  (the  roughness  length 
for  temperature)  is  required,  which  is  determined  from  the  modification  of 
Verma  [33]  and  is  expressed  as 
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1  + 


1  -  B* 


(29) 


where 

B* 

A 


y(z  -  <0  1 

y(z  -  d)  -  1 


+  2  tan'*  X 


1 


|W. 

\zo  ’ 


(30) 

(31) 
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y 


1/4 


(32a) 


1  -  15  £ 

1/2  . 

.  X  * 

1  -  15  £ 

L 

L 

< 

for  unstable  conditions, 


(32b) 


for  stable  conditions,  and 


Z*  -  0.13Z, 


(33) 


for  neutral  conditions. 

In  practice,  it  is  common  to  approximate  2^  =  0.13  Zo  at  all  times  of  the  day 
including  transition  or  neutral  periods.  Values  for  kB  *  (where  k  is  Karman’s 
constant  and  B  is  the  Stanton  number  [34])  for  different  soil  surfaces  and 
veget3tion  are  given  in  figure  4. 

Having  an  estimate  of  T,  for  dry  cloudless  conditions,  Rj.,  is  computed.  From 
the  above  equations,  the  first  estimate  of  net  radiation  can  be  computed. 

3.4  Ground  Heat  Flux 

In  Rachele  and  Tunick,  [1]  the  ground  heat  flux  is  expressed  as 
c  -G-  .(r,-r,)K.5in  { 

where 

O’  =  (r/r.  *  2)  -  r^]*r,siii{T/i2[((, .  2)  -  g} 

’^la  =  ground  temperature  during  adiabatic  conditions  (approximately  1  h 
after  sunrise) 

tn  =  time  {relative  to  midnight)  of  adiabatic  conditions 
Ks 

K„  =  - for  soil  type  and  wetness 

lyfK 

Ks  =  soil  thermal  conductivity  (W  m'*  K'*) 

K  =  soil  thermal  diffiissivity  (m^  h''). 


(34) 

(35) 


\ 
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Figure  4.  kB'*  is  a  function  of  the  roughness  Reynolds  number.  Re*  —  u*Z/v  (adapted 
from  W.  Bnitsaert,  1984,  Evaporation  in  the  Atmosphere^  D.  Reidel  Publishing  Company, 
Dordrecht). 


Equations  34  and  35  are  modifications  to  the  expressions  for  ground  heat  flux 
given  by  Angus-Leppan  and  Brunner.  [12]  Table  7  contains  values 
of  thermal  constants  newly  compiled  by  Rachele  for  this  report. 

Equation  35  implies  that  the  ground  heat  flux  is  zero  approximately  2  h  after  the 
neutral  period,  and  approximately  3  h  after  sunrise.  More  recent  experience 
confirms  that  the  morning  transition  period  may  be  approximated  by  +  1  h 
and  that  G*  can  often  be  set  to  a  constant  value  (for  example,  G*  1 1.0  -  15.0). 

The  equations  necessary  for  computing  a  second  approximation  for  H,  assuming 
that  for  clear  sides  and  dry  ground  conditions  the  latent  heat  flux  is  negligible, 
are  now  available.  Let  H  =  R.  -  G  and  complete  an  iterative  process  to 
determine  a  dry  H.  In  the  second  half  of  the  computations,  cloud  cover  on  net 
radiative  fluxes,  soil  moisture  on  ground  heat  storage,  and  evaporation  and  wind 
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Table  7.  Average  thermal  properties  of  soils,  rock,  snow,  ice,  and  water* 


I 


Material 

State 

K.  (W  m-'  K  ') 

K  (m^  hrs'*) 

Ko 

Sand 

Dry 

0.277 

9.000xl(r* 

5.0 

Sand 

Moist 

2.0 

2.520x10’ 

20.8 

Sand 

Saturated 

2.2 

2.120x10’ 

24.3 

Quartz  Sand  (Medium  Fine) 

Dry 

0.^ 

7.200x10* 

5.037 

Quartz  Sand 

8.3%  Moist 

0.586 

1.188x10’ 

8.7 

Yolo  Silt  Loam 

Dry 

0.12 

4.500x10' 

2.79" 

Yolo  Silt  Loam 

Moist 

0.44 

1.860x10’ 

5.23" 

Yolo  Silt  Loam 

Saturated 

0.848 

1.730x10’ 

10.44" 

Sandy  Loam 

Dry 

0.23 

.5.760x10^ 

4.9 

Sandy  Soil 

Saturated 

2.2 

2.660x10’ 

21.8 

Sandy  Clay 

15%  Moist 

0.925 

1.330x10’ 

12.98 

CUySoU 

Dry 

0.3 

5.490x10* 

6.28 

Clay  Soil 

Moist 

1.47 

3.960x10’ 

11.71 

Clay  Soil 

Samrated 

1.58 

1.840x10’ 

18.90 

Clay  Paature 

— 

2.76 

4.320x10’ 

21.50 

Calcareous  Soil 

43%  Water 

0.712 

6.480x10* 

13.94 

(Calcium  Carbonate) 

Soil  (Generic) 

Very  Dry 

0.167-0.345 

7.2-11.0x10* 

3.19-5.37 

Soil  (Generic) 

Wet 

1.260-3.350 

1.44-3.6xl(r* 

17.0-28.6 

Mud 

— 

0.84 

7.920x10* 

15.28 

Peat  Soil 

Dry 

0.13 

5.400xl(r* 

2.65 

Peat  Soil 

Samrated 

0.85 

4.320.10^ 

20.95 

Concrete 

0.92 

1.730x10’ 

11.33 

Granite 

o*c 

2.72 

4.570x10’ 

20.60 

Limestone 

o*c 

2.00 

2.920x10* 

18.96 

Marble 

— 

2J0 

3.490x10’ 

19.92 

Sandstone 

— 

2.60 

4.070x10’ 

20.87 

Snow 

Fresh 

0.08 

3.600x10* 

2.16 

Snow 

Old 

0.42 

1.440x10’ 

5.67 

Ice 

— 

2.24 

4.180x10’ 

17.75 

Water 

4*C 

0.57 

5.040x10* 

13.00 

'These  values  represent  coarse  averages  of  values  taken  from  Oke,  [35]  Campbell,  [14]  Geiger,  [36] 
Danard,  [16]  Van  Wjk,  [37]  Lettau  and  Davidson,  [38]  and  the  Smithsonian  Meteorological  Tables.  [39] 
'These  values  are  adjusted  to  approximate  surface  values. _ 


speed  on  sensible  heat  estimates  are  conudered.  Section  4  provides  a  step-by-step 
procedure  for  computing  the  dry  and  wet  sensible  and  ground  heat  fluxes  and  the 
net  radiation  using  equations  presrated  above.  The  latent  heat  flux,  L’E,  is 
estimated  by  the  nuxlel  as  a  residual.  Latent  heat  flux  is  expressed  as 
L’E  =  R.  -  H  -  G. 
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4.  Computational  Procedure 

This  section  presents  a  flowchart  (flgure  5)  and  step-by-step  method  description 
of  the  computational  procedure.  For  the  procedure,  the  following  are  assumed 
to  be  known  at  a  geographical  site  selected  for  computing  estimates  of  sensible 
and  latent  heat  fluxes:  (1)  longitude  and  latitude  of  the  site;  (2)  day  of  year  and 
time  of  day;  (3)  reference  level  (2  m)  values  of  temperature  (T,),  pressure  (P,), 
relative  humidity  (f,),  and  horizontal  wind^)eed  (V,);  (4)  measure  or  judgment  of 
cloud  cover  (cc)  in  tenths  and  cloud  density  on  a  scale  of  0  to  3  (3  being  most 
dense);  (S)  albedo  of  the  soil,  including  vegetative  cover;  (6)  composition  and 
wetness  of  the  soil  (dry,  moist,  saturated);  and  (7)  estimate  of  cloud  type  and 
height. 

Step  1:  Compute  the  elevation  angle  <f>  =  90°  -  Z  of  the  sun  at  the  times  of 
interest  using  Woolfs  [21]  equations  (equations  7  through  13). 

Step  2:  Compute  a  first  2q)proximation  of  the  sensible  heat  flux  for  dry 
ground  and  cloudless  skies  using  Angus-Leppan  and  Brunner  [12] 

r  C  W  sin  (36) 

(for  C  =  W  =  1  in  flgure  3) 

Step  3:  Compute  the  friction  velocity,  u*,  and  the  Obukhov  length,  L 
(equations  25  and  26a  or  26b). 

Step  4:  Compute  the  potential  temperature  scaling  constant  6*  using  equation  27 
repeated  here  as 

e*  =  u'*  kgL,  (37) 


where 

=  r,(  1+0.61  g,)  (38) 

and 

q,  =  3.8x10"’  -i-  expi 5.44x  10’  f  — - ?-  (39) 

100  ”  (273.15  rJJ 

Step  5:  Knowing  T„  $•,  and  L,  compute  the  effective  ground  temperature  T,  at 
height  Zh  using  equations  28  through  32c. 


31 


32 


Step  6:  Compute  the  downward  long-wave  radiation  for  clear  skies  (equation  31) 
using  Swinbank.  [29] 

Step  7:  From  Meyers  and  Dale  [19],  compute  R,t  for  cloudless  skies  as  in 
equations  5  through  19. 

Step  8:  Compute  sunrise  and  sunset  time  using  Woolfs  [21]  equations 
(equations  13  and  14). 

Step  9:  Compute  the  incoming  insolation  after  surface  reflection,  (1  -  a)  R,t, 
from  Paltridge  and  Platt.  [27]  Equation  20  is  repeated  here  as 

a  =  a,  +  (1  -  a,)  exp  [-if(90“  -  Z)],  (40) 


where 

k’  =  0.1 

Z  =  zenith  angle 

a,  =  albedo  at  high  solar  elevations. 

Step  10;  Compute  the  upward  long-wave  radiation  Rl,  using  Yamada,  [17]  as 
in  equation  23. 

Step  11:  Compute  a  first  approximation  for  the  net  radiation  (stq)s  1  through  10) 
as 

«„  =  (!-  «)  -  (S„  -  «„).  (41) 

Step  12;  Compute  T,  -  T^,  using  =  T„  +  (Z,  -  (Z,,  +  d)) 

where 

=  the  temperature  in  degrees  kelvin  at 
y^  =  the  dry  adiabatic  lapse  rate 

=  the  reference  height  temperature  at  the  time  1„  (relative  to  midnight) 
of  neutral  conditions. 

This  temperature  difference  will  be  used  to  compute  a  first  estimate  of  the  dry 
ground  heat  fiux,  G^. 

Step  13:  Compute  dry  value  of  G  (GJ  using  an  equation  by  Rachele  and  Tunick 
(see  equations  34  and  35  and  table  7). 

Stq>  14:  For  strictly  unstable  surface  layer  conditions,  compute  a  second 
approximation  of  using 
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Recall  that  the  sensible  heat  flux,  H,  remains  constant  for  stable  conditions  while 
H  s  0  during  transition  periods. 

Step  IS:  Repeat  Steps  2  through  14  until  convergence  is  attained,  giving  dry 
values  of  H,  R^,  and  G. 

Step  16:  For  unstable  conditions,  redefine  H  as  H  =  H,,  C  W. 

Step  17:  Determine  C  and  W  from  Angus-Leppan  and  Brunner  [12]  (figure  3) 
as  a  function  of  cloud  cover,  soil  w^ess,  and  wind  speed. 

Stq)  18:  Alternatively,  for  grassy,  damp  vegetative  ground  cover,  compute  H 
using  the  Rachele/Tunick  empiricsJly  derived  function 

H  ^  H^C  W(V;),  (43) 

where  W(Vr)  is  shown  on  figure  6  (RT  modified  Angus-Leppan,  Brunner  plot). 

Step  19:  Compute  u*,  6*,  L,  T,  using  corrected  values  of  H. 

Step  20:  If  necessary,  compute  R.^  for  cloudy  skies  using  table  4  from 
Haurwitz.  [20] 

Step  21:  If  necessary,  recompute  (1  -  a)  R,«. 

Step  22:  If  necessary,  recompute  R^t  for  cloud  cover  (equation  22). 

Step  23:  Compute  Ri.t  using  the  newly  computed  value  of  T,. 

Step  24:  Recompute 

=  (1  -  a)  R,,  -  (R^,  -  R,,).  (44) 

Stq)  25:  Recompute  T,  -  T,„. 

Stq)  26:  Recompute  a  value  for  G  using  expressions  in  equations  34  and  35  with 
the  appropriate  values  for  the  thermal  constants  outlined  in  table  7. 

Step  27:  Compute  the  latent  heat  flux  L’E  as  L’E  =  R„  -  G  -  H. 


Surface  wetness  foctor 


5.  Model  Results 

The  results  of  this  study  are  based  on  three  sets  of  measured  data  (1/2-h  averages) 
collected  at  Davis,  CA,  during  summer  1966  and  spring  1967.  [40,41,42]  The 
Davis  field  site,  a  flat,  S-ha  area,  at  17  m  above  sea  level,  is  located  about  2  km 
west  of  the  main  part  of  the  University  of  California  at  the  Davis  Campus,  24  km 
west  of  Sacramento,  and  113  km  northeast  of  San  Francisco.  Data  were  taken 
when  the  surrounding  fields  were  mostly  crop  covered  and  wdl  irrigated,  giving 
homogenous  surface  conditions  with  req)ect  to  temperature  and  moisture. 
Advection  effects  were  considered  negligible.  Profiles  of  wind,  temperature,  and 
moisture  were  measured  with  transducers  at  nine  levels  from  25  to  600  cm.  Raw 
data  were  processed  to  give  1/2-h  average  profiles. 

In  addition,  1/2-h  values  of  net  radiation  and  sensible,  latent,  and  soil  heat  fluxes 
were  available.  The  terrain  at  Davis  was  relatively  smooth  and  covered  with 
fescue  grass  (average  height  10  cm).  The  soil  was  assumed  to  be  a  moist,  silt 
loam.  Two  data  sets  were  collected  during  cloudless  sky  conditions. 

Figure  7  shows  1/2-h  average  values  of  reference  level  temperature,  relative 
humidity,  and  windspeed  for  the  3  days.  Figures  8,  9,  and  10  show  the  measured 
and  model  values  of  sensible,  latent,  and  ground  heat  fluxes  for  each  day.  Figure 
11  shows  the  measured  and  modeled  net  radiation  for  each  day.  Figure  12  shows 
the  results  for  the  numerically  approximated  optical  turbulence  structure 
parameter  (in  the  visible)  for  each  day,  using  measur^  versus  modeled  fluxes  for 
inputs. 
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Figure  7.  1/2-hr  average  values  of  reference  level  (2  m)  (a)  temperature,  (b)  relative 
humidity,  and  (c)  windspeed  (Davis,  CA  data)  6-2-66  (small  dashed  line),  7-13-66  (solid 
line),  5-M7  (large  dashed  line). 
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Figure  9.  Measured  (solid  line)  and  modeled  (dashed  line)  values  of  latent  heat  for 
Davis,  CA  data  on  (a)  6-2-66,  (b)  7-13-66,  and  (c)  5-9-67. 
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njc  (PST) 

Figure  10.  Measured  (solid  line)  and  modeled  (dashed  line)  values  of  ground  heat  flux 
for  Davis,  CA  data  on  (a)  6-2-66,  (b)  7-13-66,  and  (c)  5-9-67. 


TI»€(PST) 

(c) 

figure  12.  Optical  turbulence  structure  parameter,  C*,  for  visible  wave-lengths  for 
Davis,  CA  data  (a)  6-2-66,  (b)  7-13-66,  and  (c)  5-9-67.  Measured  fluxes  (solid)  and 
modeled  fluxes  (dashed)  as  input. 


6.  Discussion 


The  model  assumes  that  the  sensible  heat  flux  remains  constant  before  sunrise  and 
after  sunset;  therefore,  dq)artures  in  model  estimates  from  the  observed  data  at 
these  times,  as  shown  in  figures  8a  and  8b,  are  not  surprising.  A  relatively 
uniform  heat  flux  with  time  from  midnight  to  sunrise,  as  suggested  by  Angus- 
Lq)pan  [31],  is  observed.  With  regard  to  figure  8c,  it  is  clear  that  the  structure 
of  Ae  observed  data  for  the  variably  cloud  covered  sky  period  has  been 
reasonably  well  rq)resaited.  Success  relies  primarily  upon  the  s^ility  to  inteipr^ 
the  cloud  factor  effects  as  shown  in  figure  3.  The  interrelation  betweoi  the 
insolation  ratios  found  in  table  4  and  the  cloud  factors  given  in  flgure  3  must  be 
considered.  Application  of  these  corrections  for  cloud  cover  depoids  heavily 
upon  good  cloud  density  information  (and  the  more  standard  cloud  base  and  cloud 
type  input  data).  Imagine  a  partly  cloudy  sky  over  a  site.  The  rays  of  the  sun 
may  be  mostly  blocked.  However,  it  is  possible  to  observe  a  clear,  unobstructed 
line  of  sight  to  the  sun,  even  if  8/10  or  9/10  cloud  cover  were  reported.  This 
basic  scenario  can  make  radiation  and  energy  balance  modeling  extremely 
difficult,  unless  more  effort  is  directed  toward  integrating  cloud  cover  and  ground 
wetness  factors  (figure  3)  with  tte  insolation  ratio  values  of  Haurwitz  [20]  given 
in  table  4. 

Agreement  between  measured  versus  modeled  values  of  latent  heat  flux  are 
illustrated  in  figures  9a,  9b,  and  9c.  The  agreement  is  good  because  L*£  is 
computed  as  the  residual  of  the  other  estimated  energy  balance  components.  A 
strong  dependency  on  figure  11c  (net  radiation)  is  implied  by  figure  9c. 

The  best  results  of  ground  heat  flux  estimates  are  shown  in  figure  10b.  From 
each  flgure  (10a,  10b,  and  10c)  it  is  clear  that  the  expressions  applied  to  the  most 
difficult  problem  need  refinemoit,  especially  before  sunrise  and  ^ter  sunset.  The 
daytime  peak  magnitudes  are  well  rq)resented,  primarily  because  of  the  proper 
application  of  the  soil  thermal  ccmstants  given  in  table  7. 

The  most  recognizable  successes  are  illustrated  in  figures  1  la,  1  lb,  and  1  Ic.  The 
equations  chosen  for  making  estimates  of  the  net  radiative  flux  were  well  suited 
to  model  constraints  and  ultimate  application  to  electromagnetic  propagation  and 
imagery.  A  sensitivity  to  modeling  the  net  radiation  under  cloud  covered  skies 
that  relies  heavily  upon  having  good  cloud  density  approximations  was  found. 
By  using  the  information  given  in  table  4,  the  calculated  shortwave  component 
was  adjusted  to  give  better  net  radiation  estimates.  It  is  more  difficult  to  model 
radiation  exchanges  when  the  cloud  field  is  not  uniform. 
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7.  Summary  and  Conclusions 

The  purpose  of  this  study  was  to  determine  the  feasibility  of  structuring  a 
radiation/energy  balance  model  that  yields  estimates  of  sensible  and  latent  heat 
fluxes,  suitable  for  imagery  and  EM  propagation  assessments,  research,  and 
ai^lications,  when  constrained  to  require  atmospheric  measurement  of 
temperature,  pressure,  relative  humidity,  and  wind  at  a  reference  height  (about 
2  m). 

The  model  presented  satisfies  the  above  constraints;  however,  one  must  also  know 
the  day  of  the  year,  time  of  day,  longitude  and  latitude  of  the  site  of  interest, 
judgment  of  soil  type  and  moisture  (dry,  moist,  saturated),  and  cloud 
characteristics  (tenths  of  cloud  cover,  density,  cloud  type,  and  approximate 
height). 

The  model  is  a  composite  of  equations,  some  purely  physical  and  well  founded, 
others  semiempirical,  and  a  few  that  are  strictly  empiri(^.  Major  contributions 
include  equations  for  estimating  the  effect  of  wind  on  the  sensible  heat  flux;  an 
expression  for  estimating  ground  heat  flux;  a  modification  of  Verma’s  expression 
for  the  computation  of  the  roughness  length  of  temperature;  and,  most  critically, 
tying  together  the  equations  and  establishing  a  calculation  procedure. 

The  three  cases  presented  in  this  Tqx>rt  were  restricted  to  using  data  from  Davis, 
CA  to  be  consistent  with  our  previous  report.  [1]  Results  are  most  encouraging 
for  determining  qualitative  estimates  of  Ae  radiation  and  energy  balance  fluxes 
and  throughout  the  day  over  fescue  grass. 

It  is  only  fair  to  say  that  the  model  should  be  used  and  refined  for  barren  surfaces 
and  other  vegetative  cover.  A  joint  data  collection  effort  with  the  U.S. 
Department  of  Agriculture,  Agricultural  Research  Service  was  carried  out  in  May 
and  July,  1992,  in  Bushland,  Texas  [43]  to  serve  that  purpose. 
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